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Using the latest data concerning dNch/dr/ at y^sjvjv = 200 GeV, we have analyzed them 
by means of stochastic theory named the Ornstein-Uhlenbeck process with two sources. 
Moreover, we display that z,. = r;/ yjljf) ~ ?)/T?rms scaling with centrality cuts are described 
by two Gaussian distributions. 

1. Introduction In previous papers, we have shown that there is an r/-scahng 
in the charged multiphcity distributions dNch/drj or {Nch)^^dNch/dr] = dn/drj at 
■^/sjVTV = 130 GeV by PHOBOS Collaboration. ^) Those distributions have been 
explained by a stochastic approach named the Ornstein-Uhlenbeck (OU) process 
with the evolution parameter t, the frictional coefficient 7 and the variance : 



dP{y, t) 
dt 



7 



d_ 1^2^ 



Piy, t) 



(1) 



Introducing two sources at ±ymax = ^^{\/snn /it^n) at t = and P{y, 0) = 0.5[5{y-\- 
ymax)+'5(2/— ymax)]) we obtain the following distribution function for dn/drj (assuming 
y ^ rf) using the probability density P{y, t) 



dn 
dr] 



1 



V^TrV'it) 



exp 



(r] + Vn 



-7*^2 



2y2(t) 



+ exp 



2V^{t) 



, (2) 



where V^{t) = {a'^/2-/)p with p = l- e-'^'^K 

On the other hand, in Ref. 3), PHOBOS Collaboration has reported new anal- 
yses at 200 GeV. Thus, we are interested in comparisons of our results in Ref. 1) 
with theoretical analyses of data in Ref. 3), in particular Eq. (2) and its Zr = rj/rij- 
(7?rms = VW)) scaling function with Zmax = '?/max/??rms and V^{t) = V'^{t)/ril 



1 rms' 



dn 1 

d^r ~ T^Wi^y 



exp 



2V?{t) 



+ exp 



{Zr 



2V?{t) 



.(3) 



First we have to examine r/ scaling of dn/drj in a semi-phenomenological point of 
view. After analyses in terms of Eqs. (2) and (3), we examine the meaning of the 
evolution parameter t with 7, assigning a physical dimension [second] to it. Finally 
concluding remarks are given. 
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2. Semi-phenomenological analyses of data In Fig. 1, we show distributions of 
dn/dr]. As is seen in Fig. 1, the intercepts of dn/dr]\ri=o with different centrahty cuts 
are located in the following narrow interval 



dn 
dr) 



0.123-0.130. 



(4) 



ri=0 



We can observe that the ry-scaling approximately holds at y^s/vjv = 200 GeV. 

Next, we should examine a power-like law in (0.5(A^part)) "'^(c^-^ch/c^'?)|»7=0 ((-^part) 
being number of participants) which is proposed by WA98 Collaborations as 



0.5(A^part) d7] 



(5) 



??=0 



As seen in Fig. 2, it can be stressed that the power-like law holds fairly well*). Using 
estimated parameters A and a, we can express c as 



c = ^:^^^x(2.05)(Ap,rt)"-^°i 



(6) 



They are shown in Table I. The intercepts show weak decreasing behavior as the 
centrality cut increases. 
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Fig. 1. dn/dr] with different centrality cuts. 
3. Analyses of data by Eqs. (2) and (3) 

3.1 Explanation of dn/dij by Eq. (2) This time ±?7max = ±5.4 are taken. The 
estimated parameters V'^{t), p, c(Th) = {l/^/2^TV'^{t)) ■ exp [-(r/maxe~'^*)^/2^^(*)] 
and ^-re shown in Table II. The results are shown in Fig. 3. To describe the dip 
structures, the finite evolution time is necessary in our approach. 

*^ According to Ref. 4), this law suggests the following picture: The incoming particles lose their 
memory and every participants contributes a similar amount of energy to particle production. 
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PHOBOS 200 GeV 
l'*-^ ^ A{N^^)"; A = 2.05, a = 0.101, 
c.c = 0.991 
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Fig. 2. Determination of the parameters A and a. The method of linear regression is used. The 
correlation coefficient (c.c.) is 0.991 (200 GeV). From data at 130 GeV we have A = 1.79, 
a = 0.103 and c.c. = 0.993. 

Table I. Empirical examination of Eq. (6) (yijvjv = 200 GeV). 



centrality (%) 



35-45 



25-35 



15-25 



6-15 



0-6 



(A^part) 

c (Ex) 
c (Eq. (6)) 



93±5 
1230±60 
0.123±0.012 
0.122±0.009 



138±6 
1870±90 
0.124±0.012 
0.124±0.008 



200±7.5 
2750±140 
0.127±0.012 
0.127±0.008 



277±8.5 
3860±190 
0.129±0.012 
0.130±0.008 



.344.5±11 
4960±250 
0.130±0.012 
0.128±0.008 



Table II. Parameter values obtained in our analyses using Eq. (2). The evolution of P{y, ymax, t) 
in Eq. (2), i.e., dn/dr] {y « r/), is stopped at minimum values of x^- Hero, 5p = 0.004-0.006 
and Sct — 0.005-0.006. (n.d.f. means the number of degree of freedom.) Values of r^rms are 
also shown. Notice that c*(Th) = c(Th) • -^cfe'^V-^i/f"^- This ratio is needed for correction for 
\rj\ > 5.4. 



centrality (%) 


45-55 


35-45 


25-35 


15-25 


6-15 


0-6 


P 


0.855±<5p 


0.861±<5p 


0.864±(5p 


0.868±(5p 


0.873±(5p 


0.876±(5p 


v'{t) 


3.62±0.26 


3.49±0.23 


3.31±0.20 


3.09±0.17 


2.95±0.15 


2.79±0.13 




780±13 


1270±21 


1930±30 


2821±43 


3951 ±60 


5050±77 


c*(Th) 
xVn.d.f. 


0.121±<5ct 


0.123±(5ct 


0.124±<5ct 


0.125±(5ct 


0.127±5ct 


0.127±(5ct 


1.07/51 


0.91/51 


0.88/51 


1.18/51 


1.06/51 


1.46/51 




2.57 


2.54 


2.51 


2.48 


2.45 


2.42 



3.2 Comparison with other approaches First we consider a problem between 
the role of Jacobian and dip structure at 77 « 0. The authors of Refs. 5) and 6) have 
explained dNch/dr] by means of the Jacobian between the rapidity variable (y) and 
the pseudorapidity (77): The following relation is well known 

dn p dn coshr/ dn 

dr] ~ Edy ~ ^i + rn^/pl + siiiY^r^dy ' 

where dn/dy = {1/ ^/2^^V'^{t)) exp [-y'^ /2V'^{t)] . It is worthwhile to examine whether 
or not dn/drj at y^sjviv = 200 GeV can be explained by Eq. (7). As is seen in Fig. 4, 
for dn/dfj in the full phase space {\'r]\ < 5.4), it is difficult to explain the rj distribu- 
tion. On the other hand, if we restrict the central region (|?7| < 4), i.e., neglecting 
the data in 4.0 < \r]\ < 5.4, we have better description. These results are actually 
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Fig. 3. Analyses of dn/dr] with centrality cuts using Eq. (2). (See Table II.) 



utilized in Refs. 5) and 6). In other words, this fact suggests us that we have to 
consider other approaches to explain the dip structure in central region as well as 
the behavior in the fragmentation region. In our case it is the stochastic theory 
named the OU process with two sources at ij/max and at t = 0. 

3.3 Zr scaling in dn/dr] distributions The values of ryrms = VJtF) are calculated 
in Table II. The Zr = r//?7rms distributions are shown in Fig. 5(a). As seen there, 
Zr scaling in dn/dr] distributions with different centrality cuts holds at ^/sj^ = 200 
GeV. To compare Zr scaling at 200 GeV with one at 130 GeV, we show the latter 
in Fig. 5(b). It is difficult to distinguish them. This coincidence means that there is 
no change in dn/dzr as colliding energy increases, except for the region of |z,.|^2.2. 

4. Interpretation of the evolution parameter t with 7 In our present treatment 
the evolution parameter t and the frictional coefficient 7 are dimensionless. When 
we assign the meaning of second [s] to t, the frictional coefficient 7 has the dimension 
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Fig. 4. Analyses of dn/drj by means of single Gaussian and Eq. (7). (a) Data in full-77 variable are 
described by V'\t) = 5.27 ± 0.26 and m/pt = 1.13 ± 0.10. The best = 20.0/51. (b) Data in 
< 4.0 are used, i.e., 4.0 < I77I < 5.4 are neglected. V'^{t) = 7.41±0.85 and m/pt = 0.82±0.13 
arc used. The best = 4.0/37. Introduction of renormalization is necessary, due to the 
Jacobian 
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Fig. 5. Normalized distribution of dn/dzr with Zr = ??/?7rms scaling and estimated pararnotors using 
Eq. (3). (a) y/s^ = 200 GeV, p = 0.867 ± 0.002, Vr{t) = 0.506 ± 0.011, xVn-d-f- = 14.6/321. 
(b) ^/JITn = 130 GeV, p = 0.854 ± 0.002, V^{t) = 0.494 ± 0.010, x^/n.d.f. = 25.5/321. Dashed 
lines are magnitudes of error-bars. Notice that {z^) = 2max(l — p) + ^ 1.0, due to the sum 
of two Gaussian distributions. 



of [sec ^ = (1/3) X 10 fm ^]. The magnitude of the interaction region in Au-Au 
colhsion is assumed to be 10 fm. See, for example, Ref. 7). t is estimated as 

t K, 10 fm/c « 3.3 X 10"^^ sec . (8) 

The frictional coefficient and the variance are obtained in Tabic III. They are com- 
parable with values [xy^ = 0.1 — 0.08 fm~^] of Ref. 8), which have been obtained 
from the proton spectra at SPS collider. 

5. Concluding remarks cl) We have analyzed dn/drj distribution by Eqs. (2) 
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Table III. Values of 7 and at ^/s]^ = 200 GeV provided that f « 3.3 x 10"^^ sec. 



centrality (%) 


45-55 


35-45 


25-35 


15-25 


6-15 


0-6 


average 


7 [fm-^1 


0.096 


0.099 


0.100 


0.101 


0.103 


0.104 


0.101 


[fm-i] 


0.817 


0.800 


0.763 


0.720 


0.696 


0.666 


0.744 




8.51 


8.08 


7.63 


7.13 


6.76 


6.40 


7.42 



and (6) . The intercepts c's do not show remarkable energy dependence*^ . 

c2) We have shown that the dip structures in dn/dr] are hardly explained only by 
the Jacobian factor from d.n/dy in Fig. 4. What the authors of Refs. 5) and 6) 
have stressed is the explanation on the dip structures of dn/drj in the central region 
(1 77 1 < 4.0) by the Jacobian factor. This situation should be changed in analyses of 
dn/dr] in the whole pseudo-rapidity region. (See Fig. 4.) 

c3) From the evolution parameter t with 7 using the assumption for the size of 
interaction (i?(Au-|- Au) « 10 fm), we have estimated the frictional coefficient which 
is almost consistent with values at SPS energies in Ref. 8) 

c4) From dn/dz^s at y^sjviv = 200 GeV and 130 GeV, we have shown that both 
distributions are coincided with each other. If there are no labels (200 GeV and 130 
GeV) in Fig. 5, we cannot distinguish them. This coincidence means that there is 
no particular change in dn/drj between ^/snn = 200 GeV and 130 GeV. 
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*' We have observed that the intercepts c arc decreasing as the colliding energy increases: c (130 
GeV) = 0.125 0.138 c (200 GeV) = 0.123-0.130. From the semi-phenomenological analyses, i.e., 
Eq. (6), the reason is attributed to the changed ratio ^(200 GcV)^^{i30 GeV) ^ -^^^ decreased 
ratio A/'ch(130 GeV)/Afch(200 GeV) 0.82. 1.15 x 0.82 w 0.94. On the other hand, in the view of 
the stochastic approach, i.e., Eq. (2), the reason is mainly attributed to the rjmax = ln(y^ivjv/w-jv) 
as 

c (200 GeV) _ y(t)<^°°> 
c (130 GeV) ~ y(t)(i30) 



'/max ^ '/max 



21/2 (^) (200) 21/2 (^) (130) 



0.94 , 



where the suffixes mean colliding energies. 
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